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ABSTRACT 

All of the solar system gas giants produce electron cyclotron masers, driven by the solar 
wind impinging on their magnetospheres. Extrapolating to the planet orbiting r Boo, various 
authors have predicted that it may be within the detection limits of the 4-meter wavelength 
(74 MHz) system on the Very Large Array. This paper reports three epochs of observations of 
r Boo. In no epoch do we detect the planet; various means of determining the upper limit to the 
emission yield single-epoch limits ranging from 135 to 300 mJy. We develop a likelihood method 
for multi-epoch observations and use it to constrain various radiation properties of the planet. 
Assuming that the planet does radiate at our observation wavelength, its typical luminosity must 
be less than about 10 16 W, unless its radiation is highly beamed into a solid angle O -C 1 sr. 
While within the range of luminosities predicted by various authors for this planet, this value is 
lower than recent estimates which attempt to take into account the stellar wind of r Boo using 
the known properties of the star itself. Electron cyclotron maser emission from solar systems 
planets is beamed, but with characteristic solid angles of approximately 1 sr illuminated. Future 
long- wavelength instruments (e.g., the Long Wavelength Array and the Low Frequency Array) 
must be able to make typical flux density measurements on short time scales (sa 15 min.) of 
approximately 25 mJy in order to improve these constraints significantly. 

Subject headings: planetary systems — radio continuum: stars 



1. Introduction 

The star r Boo is an F6IV star located 15.6 pc 
away (jPerrvman et al.lll997fl that is orbited by a 
planet with a minimum mass of 4.14 Mj (Jovian 
masses) orbiting in a 3.3 day period (semi-major 



lengths jZarka et al.lll997 ; Farrell. Desch. fc Zar 



1999; 
2005; 



Zarka et al.ll2001t lLazio et al.ll2004t Stevens! 
Griessmeier et al l 120051 IZarkalbOOa 2007 



axis of 0.047 AU, Butler et al 



period 
Il997tl . 



By analogy to the "magnetic planets" in the 
solar system (Earth, Jupiter, Saturn, Uranus, and 
Neptune) , there have been various predictions that 
Jovian-mass extrasolar planets should also emit 
intense cyclotron maser emission at radio wave- 



Extrapolating empirical relationships based on 
the solar system planets, it is possible to make 
quantitative predictions for both the characteristic 
emission waveleng th and radio luminosity for an 
extrasolar planet ( Farrell et al.l 1999t Lazio et al 



2004 IStever3l2005l iGriessmeier et aill2005l) 

For the planet orbiting r Boo, these predic- 
tions are that its characteristic emission wave- 
length should be between about 5 and 7 meters 
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(45 and 60 MHz) and that its radio luminosity 
would result in a flux density at the Earth be- 
tween roughl y 1 and 250 mJy (lLazio et a"D 120041 
StevendlioM iGriessmeier et al.ll2005l) . The lower 
flux density estimates result from treating r Boo 
as effectively a solar twin, while the higher esti- 
mates take into account its (higher) level of stellar 
activi ty resulting from it being younger tha n the 
Sun (jStevend 120051 : IcTriessmeier et~al1 120051 ) . In 
addition, for the solar system planets, variations 
within the level of solar activity can amplify the 
cyclotron maser emission process, producing ra- 
dio luminosities (and therefore flux densities) 1-2 
orders of magnitude above the nominal level. 

Indirect evidence for extrasolar planetary mag- 
netic fields comes in the form of modulations in 
the Ca II H and K lines of the stars HD 179949 
and v And, modulations that are in phase with the 
orbital periods of their planets with th e smallest 
semi- major axes (jShkolnik et al. 2005). Though 
they moni tored r Boo, n o sim ilar modulations 



Shkolnik et al. (2005) suggest that the 



were seen. 

Ca II line modulations result from energy trans- 
port related to the relative velocity between planet 
and the stellar magnetosphere. Any such Ca II 
line modulations for r Boo would then be sup- 
pressed because the star's rotation period is com- 
parable to the pl anet's orbital period (w 3 d, 
Catala et al. 2007? ). Sugg estively, though, the po- 
larization observations of I Catala et alJ (|2007l ) do 
suggest a complex surface magnetic field topology 
for t Boo, consistent with a possible interaction 
with the planet's magnetic field. 

For the solar system planets, the cyclotron 
maser emission is fairly wideband with A/ A A ~ 2 
(Av/v ~ 1/2). The Very Large Array (VLA) is 
equipped with a 4-meter wavelength (74 MHz) 
receiving system. Images can be made with 
the VLA in its more extended configurations (A 
and B) with rms noise levels of approximately 
100 mJy beam -1 . Thus, the radio emission from 
the planet may extend to wavelengths of 4 meters, 
and the planet may be detectable with current 
instrumentation. 

This paper reports three epochs of 74 MHz ob- 
servations of r Boo with the VLA. In $5] we de- 
scribe the observations; in <J3]we analyze our ob- 
servations both from the standpoint of the indi- 
vidual epochs as well as taken collectively and we 
make suggestions about the methodology for ob- 



servations with future radio telescopes; and in Sj4] 
we present our conclusions. 

2. Observations 

Our observations were conducted on three 
epochs with the VLA in its more extended configu- 
rations. Table [T] summarizes various observational 
details. The use of the extended VLA configu- 
rations provided angular resolutions comparable 
to or better than 1'; such resolution is required 
to reduce the impact of source confusion. Dual 
circular polarization was recorded at all epochs, 
and the total bandwidth was 1.56 MHz, centered 
at 73.8 MHz. 

Post-processing of 4-m wavelength VLA data 
uses procedures similar to those at shorter wave- 
lengths, although certain details differ. The source 
Cygnus A served as both the bandpass, flux den- 
sity, and visibility phase calibrator. For two of the 
epochs (1999 June 8 and 2003 September 12), visi- 
bility phase calibration was performed in much the 
same manner as at shorter wavelengths, with the 
phases determined from Cygnus A transferred to 
the r Boo data. Several iterations of hybrid map- 
ping (imaging and self-calibration) then ensued. 
For the middle epoch, the data were acquired 
as part of the VLA Low -frequency Sky Survey 
(VLSS. ICohen et alJl2007h . For these data, small 
"postage stamp" images of bright sources from 
the NRAO VLA Sky Survey (jCondon et alJll998h 
within the field of view were made every 1 min. 
The positions of the NVSS sources are determined 
at 20 cm, a wavelength at which the ionosphere 
should present only minor perturbations. The ap- 
parent positions at 4-meter wavelength were com- 
pared with the known positions from 20 cm. Low- 
order Zernicke polynomials were used to model the 
set of position offsets and thereby infer the phase 
corrections to be applied to the data. After apply- 
ing these phase corrections, the full data set was 
imaged. 

Two significant differences for the post-processing 
were the impact of radio frequency interference 
(RFI) and the large fields of view. In order to 
combat RFI, the data were acquired with a much 
higher spectral resolution than used for imaging. 
Excision of potential RFI is performed on a per- 
baseline basis for each visibility spectrum. For 
two of the epochs (1999 June 8 and 2003 Septem- 
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Table 1 
74 MHz Observational Log 



Epoch 


Configuration 


Duration 


Beam 


Image Noise Level 






(min.) 


(") 


(mJy beam -1 ) 


1999 June 8 


D^A 


281 


30 x 29 


120 


2001 January 19 


B 


75 


80 x 80 


110 


2003 September 12 


BnA 


GG 


32 x 21 


100 



NVSS J1347+1720 



Fig. 1 — The field around t Boo at 74 MHz for 
the epoch 1999 June 8. The grey scale flux range 
is linear between 240 (2a) and 900 mJy beam -1 . 
The source in the lower left is NVSS J1347+1720 
and is used to assess the accuracy of our astrom- 
etry. 

ber 12), RFI was identified and excised manually; 
for the second epoch (2001 January 19), RFI was 
identified and excised on an automated basis by 
fitting a linear function to each 10-s visibility spec- 
trum and removing outlier channels. The large 
field of view (11° at 74 MHz) mean that the sky 
cannot be approximated as flat. In order to ap- 
proach the ther mal noise limit, we used a polyhe- 
dron algorithm (jCornwell fc Perlevlll992h in which 
the sky is approximated by many two-dimensional 
"facets." 



Figures [I] [3] present 74 MHz images of the field 
around r Boo at the three epochs. 



Fig. 2. — As for Figure Q] for the epoch 2001 Jan- 
uary 19. The grey scale flux range is linear be- 
tween 220 (2a) and 1300 mJy beam -1 . 



3. Discussion and Analysis 

We have used a number of different methods in 
order to assess the limits on the presence of any 
emission at the location of r Boo. We begin with 
an assessment of the astrometric accuracy of the 
images. 

3.1. Astrometric Accuracy 

Ionospheric phase fluctuations result in refrac- 
tive shifts in the apparent positions of sources 
within the images. In all of the images (Fig- 
ures HH3J) is the source NVSS J1347+1720. The 
position of this source is determined at a shorter 
wavelength (20 cm or 1400 MHz) at which the 
ionospheric-induced position shifts are unimpor- 
tant. For each epoch, we have fit an elliptical gaus- 
sian to NVSS J1347+1720 and determined the po- 
sition from that fit. Uncertainties, in either right 
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Fig. 3.— As for Figure Q] for the epoch 2003 
September 12. The grey scale flux range is lin- 
ear between 200 (2a) and 800 mJy beam -1 . 

ascension or declination, are between 0'.'3 and 6". 

In addition, r Boo has a measured proper mo- 
tion of ji a = —480.34 ± 0.66 mas yr _1 and fj,g = 
54.18 ± 0.40 mas yr -1 , meas ured in the epoch 



1991.25 (jPerrvman et al.lll997fh Over the decade 
between the epoch of the proper motion determi- 
nation and the measurements reported here, we 
expect t Boo to have moved by no more than 
about 5". Adding these uncertainties in quadra- 
ture, we expect that the position of r Boo in 
these images should be uncertain by no more than 
about 8". The size of the beam (point spread 
function) ranges from about 25" (epochs 1 and 3) 
to 80" (epoch 2). We conclude that our astrome- 
try is accurate to a fraction of a beam width. 

3.2. Single-Epoch Limits on the Planetary 
Radio Emission 

In no epoch do we detect statistically signifi- 
cant emission at the location of r Boo. Our first 
estimate for the upper limit to any emission from 
the planet orbiting r Boo is obtained from the 
rms noise level in each image (Table [l}. Because 
our astrometry is accurate to better than a beam 
width, we set an upper limit of 2.5<r at each epoch. 
The resulting upper limits are 250-300 mJy. 

As a second means for setting an upper limit, 
we use the brightest pixel within a beam cen- 
tered on the position of r Boo to estimate the 
flux density of any possible radio emission from its 
planet. This estimate takes into account the back- 



ground level determined by the mean brightness 
in a region surrounding the central beam. Upper 
limits determined in this manner range from 135 
to 270 mJy. 

As a final means for setting an upper limit, 
we have co-added ("stacked" or superposed epoch 
analysis) the three images, a technique has been 
used with great success to find weak sources in di- 
verse data sets (e.g., sou rces contributing to the 



hard X-ray background, Worslev et al.1 2005: 



tergalactic stars in galaxy clusters, IZibetti et al 
2005) . Experience with 74 MHz images shows that 
the rms noise level in an image produced from the 
sum of N images is VN lower, as expected if the 
noise in the images is gaussian random noise. Ob- 
viously, by co-adding images, we will be less sensi- 
tive to a rare, large enhancement in the planetary 
radio emission caused by a temporary increase in 
the stellar wind flux, but we will be more sensi- 
tive to the nominal flux density, particularly for 
the h igher estimate s for the planet's flux density 
(e.g.. lStevensll2005[) . The (la) noise level in the 
co-added image is 65 mJy beam -1 , a value consis- 
tent with that expected if the noise in all of the 
images is gaussian. We do not detect any source 
at the position of r Boo at the level of 165 mJy 

(2.5a). 

Table [H summarizes these limits. At any 
given epoch, the upper limits on radio emission 
at 74 MHz from the planet orbiting r Boo range 
from an optimistic 135 mJy to, more conserva- 
tively, 300 mJy. On average, the planet's flux 
density is not larger than 165 mJy. We convert 
these flux densities to luminosities assuming that 
the emission is broadband (37 MHz = 1/2 of the 
observing frequency) and that the planet radiates 
into a solid angle of 1 sr (see below). 

This planet does not transit its host star. While 
the planetary magnetosphere might be large, the 
typical emission altitude for an electron cyclotron 
maser for solar system planets is approximately 1— 
3 planetary radii. Thus, these limits should apply 
regardless of orbital phase. 



3.3. Likelihood Estimates for Multi-Epoch 
Planetary Radio Observations 

Earlier wo r k by Bastian. Dulk. Sz Leblanc ( 2000l ). 



[Farrell et al 
iFarrell et al 




Rva bov et al 



(and 33~2 



(2004) 
reported 



and 



4 



Table 2 

74 MHz r Boo Observational Limits 



IVictnoa 


Epoch 


(mJy) 


Limit 

(xlO 16 W) 


rms noise level 


1999 June 8 


300 


2.6 




2001 January 19 


275 


2.4 




2003 September 12 


250 


2.1 


brightest pixel 


1999 June 8 


270 


2.3 




2001 January 19 


134 


1.1 




2003 September 12 


189 


1.6 


stacked image 




165 


1.4 



Note. — uminosity limits adopt a distance of 15.6 pc and 
assume an emission bandwidth of 37 MHz (= 74 MHz/2). 



detections fro m a single observation of a p l anet. 
(In addition, lYantis. Sullivan fc Ericksonl Il977 



and IWinglee. Dulk. fe Bastianl Il986l conducted 
blind searches for extrasolar planetary radio emis- 
sion, but it is not clear that they observed any 
star now known to be orbited by a planet, and 
we believe that they observed each star for only 
a single epoch.) In addition, r Boo has been 
observed at 150 MHz wit h the Giant Metre wave 
Radio Telescope (GMRT. iMajid et allExih . the 
results from which will be reported elsewhere. 
Here we have reported non-detections from multi- 
ple observations of r Boo. Both in the context of 
these observations as well as from the standpoint 
of future observations, to what extent can mul- 
tiple non-detections of a planet be used to place 
constraints on its radio emission? 

The expected flux density fo r the radio emission 
from an extrasol ar planet is (jFarrell et al. 1999; 
Lazio et al1l2004l ) 



rad 



AMD 2 ' 



(1) 



where P rac i is the luminosity or radiated power 
from the planet, Au is the emission bandwidth of 
the radiation, O is the solid angle into which the 
radiation is beamed, and D is the distance of the 
planet (or its host star) from the Sun. In turn, 



the luminosity or radiated power and the emis- 
sion bandwidth can be related to various plane- 
tary properties (e.g., mass and rotation rate of 
the planet). The standard practice has been to 
use empirical laws from the solar system to make 
predictions for P ra d and Au. 

Planetary radio emission has a characteristic 
wavelength A c or frequency v c , which is related 
to the cyclotron frequency at the lowest altitude 



at wh ich emission is able to escape (jFarrell et al 
1999). In turn, this characteristic wavelength is 



related to the magnetic moment or the magnetic 
field strength at the planet's surface. The emission 
is typically broadband, with A C /AA = v c / Av ~ 2. 
In making predictions, it has been assumed that 
extrasolar planetary radio emission is compara- 
bly broadband, and, more importantly, that any 
observational searches would be carried out at a 
wavelength near A c (frequency near v c ). From an 
observational standpoint, we shall incorporate a 
factor to take into account the possibility that a 
search may not have been conducted at an op- 
timum wavelength. Guided by the spectrum of 
Jupiter's emission, we use 

fu(v,v c ) = < (2) 
[0, v > 2v c . 

Although the step function nature of jviy^v^ 
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may seem artificial, Jupiter's spectrum does cutoff 
sharply around 40 MHz, which is approximately 
2v c . 

Combining equations ([T]) and ©, and making 
the assumption that Av ~ v c /2, the expected flux 
density of an extrasolar planet when observed at 
a frequency v is 

2 P r d 

S(P Iad ,v c ,D,;D,i/) = ——^f v (v,v e ). (3) 

We take P ra d, v c (A c ), and il to be planetary pa- 
rameters with the objective of constraining them 
from observation. The distance D to the planet 
(or its host star) is typically known from other 
means. 

Our methodology for searching for radio emis- 
sion from an extrasolar planet has been to utilize 
a radio interferometer to make images of the field 
surrounding the planet. In the absence of a source, 
the pixels in a thermal-noise limited image from a 
radio interferometer have a zero-mean normal dis- 
tributiorO] with a variance of c 2 , so that the prob- 
ability density function (pdf) of obtaining a pixel 
with a noise intensity between TV and N + dN is 

f N{N ) = ^=e~ N2 ^ 2 . (4) 

We adopt a signal model of I = S + N for the 
intensity at the location of the planet, where N is 
the (thermal) noise in the image and S is the flux 
density contributed by the planet. We assume that 
S is constant over the duration of the observation. 
If this is not the case, then the pdf of S would have 
to be incorporated into this analysis. In practice, 
perhaps the best that one could do would be to 
use what is known about the radio emission of the 
giant planets in the solar system to develop an 
appropriate pdf. For the purposes of this analysis, 
we shall treat S as a constant, which has the effect 
that we will be placing constraints on the mean 
level of radio emission from a planet. 

A detection occurs if the pixel intensity exceeds 
some threshold, I t . Thus, 

p s (I>h\S)= 'e-C-^/^dl. (5) 

J I t (7V27T 

1 The assumption of a zero-mean distribution depends upon 
the image having been made without a so-called "zero- 
spacing" flux density, that is, without a measurement of 
the visibility function at a spatial frequency of zero wave- 
lengths. This is the case for the images we analyze here. 



For simplicity, we define x = I /a and s — S/a. 

Equation (J3j> suggests that small values of fi 
would produce large flux densities. However, small 
values of Q also imply that the radiation beam is 
unlikely to intersect our line of sight. We charac- 
terize the probability of intercept as 

Pn(«) = (6) 

Z7T 

For the solar system planets, the emission is 
beamed from both the northern and southern au- 
roral regions. We presume that the line of sight to 
a planet cannot intercept both regions simultane- 
ously, so that the the appropriate denominator for 
this probability is a solid angle of 2ir sr. Clearly 
more complicated functions are possible, but we 
do not believe them warranted at this time. The 
full probability of detection is then 

p(P rad ,ft, v c ) = PsPn, (7) 

where we have made explicit the planetary param- 
eters that we seek to determine. 

Suppose we have n observations of a planet with 
m detections, with m = describing the current 
observational state. Can we place any constraints 
on the factors in equation ((3|)? 

Consider first a scries of trials in which the 
probability of detecting the planet in any single 
trial p is the same. This case corresponds to one 
in which the observations are essentially identical. 
Then the probability of detection in the several 
trials is given by the binomial probability 

P(p;m,n)= ( n )p m (l-p) n ~ m . (8) 
\m J 

The case of current interest would be that for 
which m = 0, 

P(p-0,n) = (l-p) n . (9) 

In any actual observational case, p will likely 
vary from epoch to epoch, most likely because the 
images will have different noise levels a. This is 
certainly the case if the images are obtained from 
different instruments (e.g., the VLA vs. the Giant 
Metrewave Radio Telescope). Even images from 
the same instrument can have different noise lev- 
els, though, depending upon the prevalence of RFI 
during the observations, the number of antennas 
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used, the duration of the observation, different el- 
evations at which the source is observed, etc. ( Ta- 
ble HJ. 

We assume that the observations are indepen- 
dent, that is that the probability of detecting the 
planet in any given observation is independent of 
the other observations. This assumption is cer- 
tainly warranted from the observational stand- 
point that the noise level a is independent from ob- 
servational epoch to epoch. Then the joint prob- 
ability of detection is 



JV 



N 



V 



n*=n 



i=l 



Pi (i-Pt) 



(10) 



for which the total number of trials rn and number 
of detections mi are allowed to vary from one set 
of trials to another. We give the full expression 
here, anticipating that there may be future ob- 
servations, potentially at a range of wavelengths 
(e.g., VLA vs. GMRT vs. Long Wavelength Array 
vs. Low Frequency Array, see below). For the case 
presented here, there is a single trial at each epoch 
with no detection, m = 1, rrii = 0, and N = 3, so 



(ii) 



For a set of observed threshold intensities {xt}, 
the likelihood function is given by equation (TTT]) 




Radiated Power (W) 

Fig. 4. — The (log-)likclihood function, equa- 
tion [T2J for detecting the radio emission from the 
planet orbiting r Boo as a function of its radi- 
ated power P ra( j and the solid angle f2 into which 
the radiation is beamed, assuming that it does 
radiate at a wavelength of 4 m. The gray scale 
is linear, with white indicating regions consistent 
with our non-detections and black indicating un- 
likely regions. The rms noise levels and flux den- 
sities used are listed in Tables [1] and [21 the band- 
width of the emission is assumed to be 37 MHz 
(= 74 MHz/2). The maximum occurs in the lower 
left, near P rad ~ 10 15 W and Q ~ 0.2 sr; the 
first contour is at 99% of the peak, with subse- 
quent contours at 98%, 95%, 90%, and 67%. The 
horizontal dotted line indicates the approximate 
beam ing solid angle for Jupiter ( Zarka fc Cecconi 
20041) . 



£({z t }|P rad ,fi,0 = Y[[l -p(x > x ht \s)} (12) 

i=l 

where we have made explicit the parameter depen- 
dences entering into S. 

For the observations reported here, because we 
have observations at only one wavelength, we as- 
sume that the observing frequency is sufficiently 
close to v c that f v {y, v c ) = 1. In this case, the like- 
lihood function (equation[TJ]) reduces to a function 
of only two parameters, P ra d and tt. Also, in prac- 
tice, we compute the (base 10) logarithm of the 
likelihood. 

Figure 2] shows the (log-)likelihood function, 
which we have computed with the image noise lev- 
els listed in Table [T] and the brightest pixel lim- 
its (Table [2]). The peak likelihood occurs in the 
lower central left of the plot, near P rac i ~ 10 15 W 
and il ~ 0.2 sr. The location of this peak reflects 



our choice of limits on P ra d, as lower power lev- 
els would also clearly be consistent with our non- 
detections. 

Two comments on the likelihood function are 
warranted. First, the shape of the allowed region 
reflects two competing effects. From equation ([3]), 
the quantities P ra d and Q are degenerate. The 
planet could radiate intensely but be beamed into 
a narrow solid angle, with a low probability of de- 
tection, or it could have a wide beaming angle but 
with only modest luminosity. The competing ef- 
fect is that, as the beaming angle becomes smaller, 
the probability of it intersecting our line of sight 
becomes progressively smaller. 

Second, the upper limit on P ra d in Figure [4] is 
motivated by pre dictions. A cr ucial element of 
the prediction by IStevensI ( 20051 ) is that the stel- 
lar mass loss rate was predicted on the basis of the 
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star's X-ray flux. It now seems that stellar mass 
loss rates change character at a char acteristic X- 



ray fl ux of about 8 x 10 5 erg cm s 1 (jWood et al 



2005). Below this value, stellar mass loss rates are 
highly correlated with X-ray flux; above this value, 
the mass loss rates drop by an order of magnitude 
and remain relatively constant with increasing X- 
ray fluxes. The X-ray flux of r Boo is essentially 
equal to the value at which stellar mass loss rates 
change character. Thus, its mass loss rate could 
be ov erestimated by as much as an order of magni- 
tude. Stevens! ( 2005 ) predicts that the planet's cy- 
clotron maser strength should scale with the mass 
loss rate as M 2 / 3 . Thus, even if it is the case that 
the mass loss rate was overestimated by a factor 
of 10, the radiated power from the planet could 
still be a factor of 5 above that which would be 
estimated assuming a solar mass loss rate. 

For our observations, assuming that the planet 
does radiate at a wavelength of 4 m, wc conclude 
that the upper end of the predicted power levels is 
allowed only if the beaming angle is <gC 1 sr. The 
cyclotron maser emission from the solar system 
planets is characterized by solid angles f2 > 1 sr, 
typically emitted in a hollow cone, with a fairly 
wide opening angle (~ 90°) and a finite thick- 
ness. Recent analysis of Cassini observations of 
Jupiter during the spacecraft's cruise phase in- 
dicate that its decameter radiation illuminates a 
cone of half-width opening angle of 75° and a 
thickness of 1 5°, with an equivalent solid angle 
£1 ~ 1.6 sr (jZarka &; Cecconil 120041) . In some 
cases, for example the Earth's auroral kilomet- 
ric radiation, the beaming angle can be f2 > tt sr 
(jGreen fc Gallagherlll985l ). 



For beaming angles comparable to those seen 
in the solar system (fi ~ 1 sr), the planet must 
radiate less than about 10 16 W or we would have 
a reasonable probability of detecting it. Luminosi- 
ties P rfl ri < 10 16 W are within the predicted range 



(lLazio et al.ll2004 ; Stevensfl2005 : Griessmeier et al 



2005), though below the most recent estimates 
that attempt to incorporate what is known about 
the star's stellar wind strength. For compari- 
son, Jupiter's nominal radio luminosity is of order 

10 10 W. 

Bast ian et al. (|200d) also subdivided each ob- 



servation of a single star into scans as short as 10 s 
(the shortest nominal integration time provided by 
the VLA). The motivation is that Jupiter's emis- 



sion is quite variable, so, if the extrasolar planet 
radio emission is similar, subdividing the observa- 
tion offers the possibility of detecting a "burst," 
which might be quite short in duration and "di- 
luted" if the ent ire (> 1 hr) observa tion is con- 
sidered. Though Bastian et al. ( 2000h did not ob- 
serve t Boo at 4-m wavelength, we have considered 
the impact of subdividing the observation. The 
obvious benefit of subdividing is that the effective 
number of observations increases; the disadvan- 
tage is that the noise level does as well. 

We assume similar observing parameters as ob- 
tained in our observing programs, namely a point- 
source noise level a ~ 0.1 Jy obtained in 1 hr for 
observations toward a star like r Boo. We find 
that sub-dividing the observation into 10 scans 
(e.g., each of 5 min. duration), which increases 
the noise level to a w 0.3 Jy, leads to essentially 
no improvement in the constraints that can be set. 

As a final comment on our likelihood method, 
we anticipate future observations, either with ex- 
isting or future instruments (see below) might be 
at different wavelengths. A significant assump- 
tion in our likelihood analysis is that searches are 
conducted at a common wavelength and that the 
planet emits at that wavelength. In the case of 
Jupiter, for instance, its electron cyclotron emis- 
sion cuts off sharply shortward of approximately 
7.5 meters (y > 40 MHz). Presuming that similar 
processes operate in the magnetospheres of extra- 
solar giant planets, observations at different wave- 
lengths may not be equally constraining. Simi- 
larly, one might wish to take into account what is 
known about the beaming angles from the solar 
system planets. It would be relatively straight- 
forward to incorporate this prior information and 
extend our likelihood analysis to a Bayesian for- 
mulation. In that case, appropriate priors would 
have to be specified for the emission wavelength, 
beaming angle, and radiated power. 

3.4. Future Observations 

There are a number of next-generation, long 
wavelength radio instruments under development. 
Notable among these are the Low Frequency Array 
(LOFAR) and the Long Wavelength Array (LWA). 
If deployed as intended, both promise to provide 
sensitivities at least an order of magnitude larger 
than the 4-meter wavelength VLA at comparable 
wavelengths (A > 3 m or v < 100 MHz). Here 
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we consider whether additional observations with 
current instrumentation is preferable to the oper- 
ation of these future facilities. 

We simulated two sets of flux density measure- 
ments toward a star like r Boo. For both mea- 
surements the typical flux density measurement 
was taken to be a gaussian random variable with 
a specified mean and variance and the rms noise 
level was taken to be a factor of 2.5 times smaller. 
The first set of measurements had noise levels and 
flux density limits similar to those measured here, 
but we considered the impact of having 15 mea- 
surements rather than just 3. Not detecting plan- 
etary radio emission in such a data set would not 
improve significantly on the radiated power con- 
straints (at most by a factor of a few) but would 
place increasingly severe constraints on the beam- 
ing angle at large power levels. 

We have also considered the requirements on 
an instrument in order to improve the constraints 
significantly. Figure [5] illustrates the constraints 
placed by 5 measurements with an instrument ca- 
pable of obtaining flux density measurements]! of 
approximately 25 mJy with rms noise levels of 
approximately 7 mJy beam -1 . Obtaining signif- 
icant (an order of magnitude or better) improve- 
ments on the powers radiated by extrasolar plane- 
tary electron cyclotron masers will require mea- 
surements at this level. One significant advan- 
tage that future long-wavelength instruments are 
likely to have over current instruments, though, is 
that future instruments are likely to have a multi- 
beaming capability. Consequently, there is likely 
to be significantly more time for observation and 
a much larger number of measurements (or even 
detections!) will be obtained. 

4. Conclusions 

We report three (3) epochs of 4-meter wave- 
length (74 MHz) observations of t Boo with the 
VLA. Our objective was to detect the electron 
cyclotron maser emission from the planet orbit- 
ing this star, assuming that the planet does pro- 
duce the equivalent of Jovian radio emissions. In 




2 As in lLazio et al. we assume relatively short obser- 

vation durations. Given the bursty characteristic of solar 
system planetary radio emission, it is not clear that the 
typical radio astronomical practice of conducting long in- 
tegrations is appropriate. 



1Q 15 1Q 16 

Radiated Power (W) 

Fig. 5. — As for Figure [H but for a next gener- 
ation, long- wavelength instrument (e.g., the LWA 
or LOFAR) capable of flux density measurements 
in the range of 25 mJy with an rms noise level 
near 7 mJy beam -1 . In this case, five measure- 
ments were simulated. 



none of our 3 epochs did we detect emission at 
the location of r Boo. For a single epoch (Ta- 
ble [2]), our limits on its emission range from 150 
to 300 mJy, equivalent to a range of luminosities 
of approximately 1-2.5 x 10 16 W. We have also 
stacked ( "co-added" ) the images to produce an av- 
erage limit of 165 mJy, equivalent to a luminosity 
of 1.4 x 10 16 W. 

We have developed a likelihood method to con- 
sider multi-epoch measurements. Our likelihood 
function depends upon three planetary parame- 
ters, the planet's radiated power level P ra d, the 
beaming solid angle of the electron cyclotron 
maser f2, and the characteristic wavelength of 
emission A c . Under the assumption that the planet 
does radiate at our observed wavelength, the pa- 
rameters P ra( j and £1 are degenerate (Figure 01 . 
We find that the typical radiated power must be 
less than about 10 16 W unless the beaming solid 
angle is f2 <C 1 sr, which would be considerably 
smaller than the typical value for solar system 
planets. Power levels Prad < 10 16 W a r e within 
the predicted range flLazio et al. 12004 IStevensI 
2005c iGriessmeier et al.l I2005I ). being comparable 
to or below the most recent estimates that at- 
tempt to incorporate what is known about the 
star's stellar wind strength. 

We have also considered future long-wavelength 
instruments, such as the LWA and LOFAR. We 
find that to improve significantly upon our con- 
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straints, future instruments will have to be able 
to measure typical flux densities of approximately 
25 mjy. 

We thank P. Zarka for illuminating discussions 
about solar system planetary emissions, J. Cordes 
for helpful discussions about likelihood functions, 
students in the Thomas Jefferson High School for 
Science and Technology astronomy lab for assis- 
tance with assessing the utility of stacking, and 
the referee for several helpful comments about the 
presentation and content. The National Radio 
Astronomy Observatory is a facility of the Na- 
tional Science Foundation operated under coop- 
erative agreement by Associated Universities, Inc. 
Basic research in radio astronomy at the NRL is 
supported by NRL 6.1 Base funding. 
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